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Abstract 

G\ ; 

"*^h , We consider the behaviour of open quantum systems in dependence on the coupling to 

one decay channel by introducing the coupling parameter a being proportional to the average 
degree of overlapping. Under critical conditions, a reorganization of the spectrum takes place 
which creates a bifurcation of the time scales with respect to the lifetimes of the resonance 
states. We derive analytically the conditions under which the reorganization process can 
be understood as a second-order phase transition and illustrate our results by numerical 
investigations. The conditions are fulfilled e.g. for a picket fence with equal coupling of the 
states to the continuum. Energy dependencies within the system are included. We consider 
also the generic case of an unfolded Gaussian Orthogonal Ensemble. In all these cases, the 
reorganization of the spectrum occurs at the critical value a cr it of the control parameter 
globally over the whole energy range of the spectrum. All states act cooperatively. 
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1 Introduction 



Recently, the properties of open quantum systems have been studied with a renewed interest 
in the framework of different approaches. Mostly discussed is the restructuring of the systems 
taking place at high level density under critical conditions and the resulting formation of dif- 
ferent time scales in terms of lifetimes of resonance states. The reorganization occurs if the 
degree f /D of overlapping reaches a critical value (f is the average width obtained by averaging 
over all M resonance states in a certain energy region and D is the mean level distance). It is 
investigated for resonance phenomena in nuclei ||, [2], [3|, [|, ||, |6| , atoms [0, ||] and molecules |J . 
In the meantime it has been considered also in other systems such as e.g. quantum dots (T(J 
and microwave billards |ll|] . The number M of resonance states is usually much larger than the 
number K of open decay channels. 



A powerful method in describing the properties of an open quantum system is the projec- 
tion operator technique. It allows us to investigate, in a direct manner, the corrections to the 
many-particle states in a subspace of the full Hilbert space which arise from the coupling to 
the orthogonal subspace. Originally the idea of using the projection operator technique arose in 
nuclear physics. The influence of discrete states as well as closed decay channels (Q subspace) 
onto the wavefunctions of the open decay channels (P subspace) is considered in the Feshbach 
unified theory of nuclear reactions [12]. This phenomenological model describes well the proper- 
ties of nuclei at an excitation energy of about 10 MeV by employing statistical assumptions on 
the states of the Q subspace. In contrast to this, the ground and low-lying states of nuclei are 
described well by the wavefunctions of the discrete states (Q subspace), considering the influence 
of the decay channels (P subspace) in an approximated manner. 



The properties of an open quantum system are described well by using the latter division of 
the whole function space: the Q subspace contains the discrete states of the system while the 
P subspace consists of open as well as closed decay channels In studying the restructuring, 
we are interested in the properties of the states of the Q subspace modified by their coupling to 
the P subspace playing the role of an environment. This method can be used for a wide class 



of open quantum systems [13]. 



The question whether this restructuring may be considered as a phase transition of second 
order is put in ||] but not considered in detail up to now. A possible analogy to the formation 
of laser light is investigated numerically in 0]. As in the case of the laser, a control parameter 
a can be defined which is proportional to T/D. The information entropy changes rapidly in a 
relatively small region of the control parameter in both the laser ]l4|] and the open quantum 
system ||. 

In other investigations || it was realized that the avoided crossing of two neighbouring reso- 
nance states which are coupled to one common channel is the basic process of the restructuring 
observed globally in the system. As soon as two resonances start to overlap each other, their 
interaction via the continuum can no longer be neglected. As a function of the coupling to a 
certain decay channel the two resonances approach each other in energy up to a certain min- 
imum distance in the complex energy plane at a = a^rit- The avoided crossing is reflected in 
the wavefunctions of the two resonance states. The biorthogonality reaches its maximum if 
a — > a cr it , and vanishes if a — > and a —> oo || . As a function of further increasing a > a cr it , 
the width of one of the two resonance states decreases (resonance trapping) while the width of 
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the other one increases further. 

The local resonance trapping can explain, indeed, the global restructuring of the quantum 
system under critical conditions. It determines also, as will be shown in this paper, whether 
the restructuring of the system takes place collectively with the simultaneous participation of 
all basis states or successively by individual trapping of resonance states. 

In the following, we will investigate this question in detail. In sect. 2, we write down the 
basic equations used in the paper. The model is formulated and the characteristic polynomial is 
given. The Hamiltonian is non-hermitian and its eigenfunctions are, generally, bi-orthogonal. In 
section 3, the properties of a system with picket-fence distributed levels coupled with the same 
strength to one decay channel is investigated in detail. The study is performed analytically for 
the limiting case of an infinite number of states as well as for a finite number. The results are 
illustrated by numerical calculations. Finally, the results obtained are discussed and identified 
with characteristic features of a second-order phase transition. The process of formation of a col- 
lective state aligned with the decay channel is discussed in detail. Its wavefunction is coherently 
mixed in the wavefunctions of all basis states including those states which are not overlapped 
by it. 

The results obtained in section 3 are underpinned in section 4 by considering some other 
level and coupling-strength distributions being more realistic than those in section 3. The 
study is performed both analytically and numerically. General conditions for the appearance 
of a second-order phase transition are formulated analytically and illustrated by the results of 
numerical calculations. As a special case, the sharpness of a phase transition is shown to be 
distorted by an imaginary part in the coupling term. The results are summarized and discussed 
in the last section. 



2 Basic equations 

2.1 The Hamiltonian of an open quantum system 

Let us consider the Hamilton operator 

H = H + V (1) 

of a many-particle system where Hq describes the mean-field, i.e. the motion of the particles in a 
finite depth potential and V is the operator of the two-particle residual interaction. A convenient 
method to solve the Schrodinger equation {H — E)^f = in the full Hilbert space of discrete 



and continuous states is to use the projection operator technique introduced by Feshbach |12f| . 
Here the whole function space is divided into two subspaces by using the projection operators 

M 

Q = J2\®k)($k\ 



k=l 

coo 



roc 

P = Y, dE%(E'))(UE') 

r J E c 



with P + Q = 1, where Q contains the many-particle discrete states | <]?/%) being solutions of 
(QHQ — Ek) \$k) = 0, and P the many-particle scattering states \(; C (E)) being solutions of 
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(PHP — E) \£ C (E)) = 0. The total Hamiltonian acting on the full Hilbertspace is split into four 
terms: H = QHQ + QHP + PHQ + PHP. 

We are interested in the properties of the Hamiltonian of the open quantum system, which 
acts on the Q subspace and carries the influence of the P space. The derivation of this Hamilton 
operator can be found in ||, |i~3|| , 

H$ Q {E) = QHQ + QHP ■ 4 +) (E) ■ PHQ . (3) 

It consists of two terms. The first one (QHQ) describes the behaviour of the closed system of 
discrete states which includes the configurational mixing due to the two body residual inter- 
action, but does not take into account the coupling to the decay channels. The second term 
gives the correction due to the coupling of the two subspaces and contains the propagator in the 
P-subspace G { p ] (E) = P[E + irj- PHP}- 1 P. 

Due to this propagator the effective Hamiltonian is energy dependent and non-hermitian. 
Its complex eigenvalues Xk(E) = £k(E) — i/2 Tk(E) give the poles of the resonance part of the 
scattering matrix 

ares • \ 

*~ x h E ~ ^k(E) (4) 

where jk c (E) = \j\p2/n (S, c (E)\V\^k) is the transition matrix element between a bound and a 
scattering state. Thereby, the complex eigenvalues get a concrete physical interpretation as 
the energy positions £jF es = Sk(£l° s ) and total decay widths T r £ s = rfc(££ es ) of a resonance state 
[||, [l3|]. The £k differ usually from the corresponding eigenvalues E^ of TL° = QHQ, i.e. from the 
energies of the states of the unperturbed system. So the external coupling to the decay channels 
causes not only the finite lifetime of the states but generally also an energy shift. 

In the following, we will restrict ourselves to an energy region in which the energy dependence 
of the Hamiltonian is small in spite of a large number M of states lying in it. Further, we consider 
a small number K of decay channels which are all open and not coupled among themselves. Then, 
the effective Hamiltonian (||) in the Q subspace is, to a good approximation, 

n = H° - iaVV + (5) 

where VV + is a hermitian operator if we consider time reversal invariance. As in Eq. (|3|) the 
first term 7i° describes the internal structure of the unperturbed system in the Q subspace. The 
second term iaVV + follows from QHP • (E) • PHQ and describes the coupling between the 
two subspaces. The parameter a, assumed mostly to be real, characterizes the mean coupling 
strength between discrete and continuous states. 



The Hamiltonian (J5j) is used successfully for the description of resonance states in nuclei |1J 
and molecules [p|. Nowadays, it is applied also to the description of resonance phenomena in 
other systems such as e.g. quantum dots [ 10 1 and microwave billards Q. 



The rank of 7i° is equal to the number M of states considered. Its non-diagonal matrix 
elements describe the configurational mixing of the discrete states. The coupling matrix V is 
a K x M matrix if the number of open decay channels is equal to K. The element Vf of V 
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describes the coupling of the discrete state i to the channel c; i = 1, M; c = 1, K. Thus, 
the rank of VV + is K. 

As long as a is small, the second term of the Hamiltonian TL can be considered as a small 
perturbation of TL°. This condition is always fulfilled if the average width f is much smaller than 
the average distance D between neighbouring resonance states. In this case, the non-diagonal 
matrix elements of TL are small, the individual resonances are isolated. Their positions and 
widths obtained from the eigenvalues Aj of TL differ only slightly from the real and imaginary 
parts, respectively, of the diagonal matrix elements of TL. 

In the opposite case of large a, the matrix VV + determines the behaviour of the system. 
Then, the rank of TL is given by K. That means, M — K states are almost decoupled from the 
continuum of decay channels and become long lived (trapped) while K states take almost the 
whole coupling strength: £^ =1 I\/2 « Q{Tr(TL)} and J2iLx+i^i ~ 0- Therefore, two different 
time scales arise at large a, see e.g. ]l[ |], ||, [J. 

Thus, a reorganization in the open quantum system takes place in the transition from small 
coupling parameters a to large ones when M 3> K. In the following, we will investigate the 
question whether and under which conditions the reorganisation of the open quantum system 
can be understood as a phase transition in the limit M — > oo. We restrict ourselves to the case 
with one open decay channel {K = 1). 



2.2 The characteristic polynomial 

We consider a system with M = 2N + 1 states coupled to one common decay channel (K = 1). 
The unperturbed eigenvalues of Tt° are denoted by Ek, k G {— N, N}, so that Ej < Ey. if 
j < k (without degeneration). The centre of the spectrum is assumed to be at Eq = without 
loss of generality. The coupling vector will be denoted by V = (w_jv, t>_i, vq, v i, ... 3 ujv). 



Due to K = 1, all column and row vectors, respectively, of VV + are linearly dependent. 
Substracting times the row from the row k, one gets the following expression for the 
characteristic polynomial, 



Pn(X) = 




which can be written as 



E-n - 




A 





E-N+l 



-iav- N+ i 



E-x - A 
— iav-x 



Aw-w+i 



\V-i 

-iavo - 
Xvx 







A —iavx 
Ex - A 





-iavM 




En — A 



N N N 

Pn{\)= II (E k -X)-ia- M 2 - II (^'-A) = 

k=-N k=-N j=-N,j^k 



= 



(6) 



Eq. (0) can be proven by induction. 



5 



According to Eq. (||), Pat (A) is the sum of two polynomials, 

P N (X) = Q N (X)-ia-R N (X), (7) 

where Qn is of the order 2N + 1 and Rn of the order 2N. If \vk\ 2 = 1 Vfc, the Qm and Rn are 
related in a simple manner, 

Rn = ~4^Qn ■ (8) 
dX 

In the limit a = 0, we find Xk = Ek = Ek V/c, i.e. the eigenvalues of TL are equal to those of TL° 
(according to the definition of the parameter a). 

The limit of large coupling strength (a — > oo) can be obtained when we rewrite the charac- 
teristic polynomial (0) as 



N 

P N (X) = ia J] (E k -X) 

k=-N 



1 N 2 1 



(9) 



The first factor of the product term is zero only at the unperturbed eigenvalues Ek of TL°. 
Therefore, for a / the solutions of Eq. (SJ) must be given by the zeros of the second factor, 
i.e. by the solutions of J- = J2 \ v k\ 2 In the limit a — > oo, there are 2N solutions lying 
at real energies: Xk £ (Ek,Ek-i) if k > and Xk £ (Ek,Ek+i) if k < where Ek is eigenvalue 
of the unperturbed Hamiltonian TL°. In the case of the picket fence with Ek = k and equal 
coupling, Xk approaches k± 1/2. Furthermore, we have exactly one complex solution at Eq = 
and ro — > oo for a — * oo. 

Let us now discuss the behaviour of the system as a function of increasing coupling strength 
a. /,From Eq. (^), we get 

^Q' N (X) -iR N {X) -ia^R' N (X) = (10) 
da da 

for the solutions of Pjv(A) = and further the differential equation 

dX iR N W 



(11) 



da Q' N (X) — iaR' N (X) 
with the initial condition Xk(a = 0) = Ek- 

For small a, Eq. (|ll|) reads 

dXk ii?Ar(Afc) 2 

da Q N (Xk) 

That means, the imaginary part of eigenvalue Afc of TL increases, with increasing a, proportional 
to \vk\ 2 while the real part of it remains unchanged, as long as a is small. 

For large a, we have 2N solutions whose imaginary part is small while the real part Ek is 
determined by Ek < £k < ^fc-i if k > and Ek < £k < ^fc+l if fc < 0, since Xk(a — > oo) = Ek- 



6 



The relevant part of i?jv(A) for the solutions of Pat (A) = is therefore T/v(A) = Ili=i(£fc ~~ ty- 
Inserting 

\ k (a)=£ k -i^ + 0(a- 2 ) . (13) 
a 

into 

= Q JV (A)-iaT JV (A) (14) 
leads in the two lowest orders in 1/a to 

N 2N 

° = n ( e j + *— - ^) - »« n^- - f * + w«) • ( is ) 

j=-N ° 1 



The solution is 

UjL-N( E j -£k) . 

Eq. (15) shows that for large coupling strengths, the decay widths of 2N states decrease as - 



with increasing a. This decrease is called resonance trapping. 

Besides these 2N solutions for large a, we have a solution at Eq = and To — > oo in the 
limit N — > oo. 

In sections 3 and 4, we will study in detail the properties of the characteristic polynomial 
rtq) by means of special cases. 



2.3 The eigenfunctions of a non-hermitian Hamilton operator 

Another value characterizing the reorganization which takes place in the open quantum system 
under critical conditions, is the mixing of the wavefunctions of the resonance states [jl], ||]. The 
mixing caused by the coupling of all the states to the common decay channels is related, in a 
natural manner, to the basic set of wavefunctions of the closed system, 

M 

'I'- (17) 

where <3?j are eigenfunctions of TL and <&!■ are eigenfunctions of TC°. The eigenfunctions <&j of the 
non-hermitian Hamiltonian TC are bi-orthogonal. The right and left eigenfunctions are defined 
by 

(W-AiM) = o 

(*i|(H-Ai) = (18) 



with the normalization 



(^m) = s i>j , (<i>;' <K) / (19) 



In our case |<3?') = (|$^)) T [||. In the following we will drop the indices r and I considering only 
the right eigenfunctions. Then the second relation of Eq. ( |l9| ) reads 



($i|*i> > 1 (20) 
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and i ^ j, is a complex number, generally. 

A good numerical measure for the strength of mixing is the number Nf of principal compo- 
nents in the eigenfunction <I>j. For its definition we are using the quantity 

h S = ■ (21) 

2^1=1 \ a tl\ 



Then, the number of principal components can be calculated as 

1 

W 



^-SeEHv*- (22> 



The value of Nf can be understood as a measure of (external) collectivity of the resonance state 
In the limiting case of equal mixing of the state i with all states j, bij = \j\f~M V j, we 
get Nf = 1 (maximum external collectivity). In the opposite case (no external collectivity) we 
have fey = 6 itj and Nf = 1/M. Generally, 1/M < Nf < 1. 



Further, we introduce the value 

M 

M 



i=l 



which characterizes the degree of non-Hermiticity of TL according to Eq. (po[). It is a function 
of a and B = 1 if Ti is hermitian. 



3 The ideal picket-fence distribution 

Let us consider first the simple case of a picket-fence distribution of M = 2N + 1 levels which 
are all coupled with the same strength (" ideal picket-fence distribution" ) to the continuum con- 
sisting of one decay channel (K = 1). The advantage of this simple model is that analytical 
studies can be performed. 



3.1 Analytical study for the limiting case TV — ► oo 

Suppose Ek = k und \vk\ = 1 V/c. Then Eq. ([]) reads 

Pn(X) = Q N (X)-ia-R N (X) 

N N N 

= n (k-x)-ia- x: n ^- x ) ^ 

k=-N k=-N j=-N,j^k 

and the relation (|8|) holds. In order to consider the limit N — ► oo, we divide Q n by a convergence 
ensuring factor, 

jv-,oo _J]fc =1 -(A;) 2 jg jg ii fc vr 
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Then the characteristic polynomial reads 

P(A) = sin(vrA) + ma cos(vrA) = . (26) 



Denoting the complex eigenvalue of W by A = £ — and splitting Eq. ( p6[ ) into its real and 
imaginary part we get (for real a): 



cos(7r£) e 7rr (l — 7ra) — (1 + air) 







sin(7r£) e nL (1 - vra) + (1 + avr) =0 (27) 

Since the two functions cos(x) and sin(x) never vanish for the same argument x, we have to 
consider two different cases: 



1. sin(-7r£) = 0=^£ = neZ and 



e^ r = 1 + 12 (28) 
1 — 7ra 



has a real solution T for a < - only. For small a, we have therefore 

r = Ii„(l±^) ,29) 

and r — ► — — In e for a = ~(1 — s) and e — > 0. 
2. cos(7TiS) = 0=>£ = n + iforraeZ and 

_ r 7ra + 1 

7TO — 1 

The last equation can be fulfilled only for a > ^. For large a it is therefore 



(30) 



r = Im(^) (3i) 

and r —y — - In e for a = -(1 + e) and e — > 0. 

As a result: The widths of all the states increase up to infinity as a function of increasing a. 
The singularity at the critical point a cr it is determined by ln(e). It is logarithmic. 

Further, the energetical positions of the states remain unchanged at the unperturbed ener- 
gies Ek = k of the system (eigenvalues of TLq) up to a — > —. At a cr it = ^, the real part £k of 
2A^ eigenvalues (all k 0) of 7i jumps from k to k — ^ if k > and from k to k + ^ if k < 0, 
respectively. As a function of further increasing a, the imaginary part of the eigenvalues of the 
2A^ resonance states (all k but k = 0) decreases first as ln(o;) while it approaches zero as — for 
a — > oo according to Eq. (||). 

In order to study the behaviour of the state in the centre of the spectrum at the energy 
So = we consider only the highest-order terms of A in Eq. (0): 

X 2N+1 + ia(2N + l)X 2N = . (32) 
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For large a (a 3> ~), the state corresponding to the solution A = —ia(2N + 1) lies at £ = 
and its width increases linear with a. 

Summarizing the results, we state the following: In spite of the fact, that the coupling param- 
eter a enters the equations (||) and (]?]) linearly, the imaginary parts of the complex eigenvalues 
show a singularity at the finite value a = — . For larger couplings, a clear separation of the time 
scales with respect to the decay widths of the resonance states occurs. This happens also in the 
case of an infinitely extended spectrum. This is not a local effect of a locally broad resonance in 
a restricted energy region, but it is produced by the whole system in a collective manner. All 
basic states, independent of their energy position, act cooperatively. 



3.2 Widths at the critical point for finite N: analytical study 

The sum of the widths of all states is, in our simple example with equal coupling strengths, 
given by 



Q{Tr(n)} = J2Y = ^(2N + l) 



(33) 



It is Tr(7i) = const (a). Thus, Q{Tr(7i)} should be a smooth function of a not only far from 
the critical point but also near to it in spite of the divergence of the widths for N — ► oo at a = — 
(see subsection 3.1). In the following, we will proof this statement. 



First, let us consider the eigenvalues of TL for finite N . In this case, we have 



N 



Pn(X)= II (fe-A) 

k=~N 



N 



1 — ia 



j=-N ■ 



(34) 



instead of the simple Eq. ( |26|) holding for N —* oo. As discussed in relation with Eq. @, the 
solutions of Pjv(A) = follow from 1 — ia YljL-N J^a = ^' Here, we are interested in the differ- 
ence between the solutions obtained for finite N and those for iV 



oc 



It holds 



N 



= l-ia 



1 



k- A 



k=-N 

where the correction term is given by 

1 



2A £ 



k=N+l 



k 2 - A 2 



1 + iair cot(7rA) + 2iaA } 



1 



k=N+l 



k 2 - A 2 



2A 
In 



dx 



N+l/2 X" 

y + i. 



A 2 



In 



dy 



N y 2 -1 



1-1/2/ 

1 + Vv 



ln(l _ _ + (y 2 )) 

y 



■ oo 


" 2" 


oo 


\ 

A 


. y. 


A' 

A 



2A 

77 



(35) 



(36) 



under the assumption 1/y = X/N <C 1. This condition is fulfilled, to a good approximation, in 
the centre of the spectrum. Splitting Eq. ( |35"1) into its real and imaginary parts (with A = £ — |r), 
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one arrives at 



= 1 — Q7T 



sinh(-7iT) 



a- 



cosh(TiT) - cos(27r£) N 



for the real part. Here the identity 



cot(x + iy) 



sinh(2x) — isin(2y) 



cosh(2y) — cos(2x) 
is used. The equation for the imaginary part (£ / 0) reads 







an 



sin(27r£) 



+ 



2a£ 



cosh(vrr) - cos(27r£) N 



from which we get 



2£ 2£ 

= sin(27r£) H cosh(7rr) cos(27r£). 

ttN ttN 



An estimation for the upper limit of the widths T of the states at a = - leads to 

Nit 



7rr = arcosh 



cos(2-7rc ) — — — sin(27rc ) 

Z6 



In 



\£\ 



Here, we have used N/\£\ 3> 1 which is fulfilled only in the centre of the spectrum. Thus, 

r. N i Nn 

2 (a = acrit) " 27 ln W 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



which holds for every one of the 2N states (for all k but k = 0) at the critical point. It means, 
r < In N for a -► ^ for all N. 



Using Eq. (p2[) , one gets the following estimation for the trace of the imaginary part of H' t 



off 



at a 



i . 



N 



j=-N 



N 



1 , NTT 



V — « 2 / dE — ln— = -N ]n(x)dx = -N [xla(x) - x\l 
— ' '" 2tt E " 



o 



1/tt 



,1/tt 







N, , 1. N , . . . 2N 2N + 1 
— (1 - In -) = —(1 + ln(vr)) » « — 

7T 7T 7T 7T 7T 



(43) 



The comparison of Eqs. (^) and (43) shows that Eq. ( ^3| ) holds also at the critical point. 
This means, the singularity of the decay widths V at the critical point occurs such that the sum 
rule J2i Fi = const (a) is fullfilled also for a — ► a cr ;t and iV — > oo. 



in leading order in iV by integrating Eqs. (41) over the interval (—1/2, 1/2): 



At the critical point, the width Tq of the state in the centre of the spectrum can be estimated 

(44) 



r , 

— (a = a 
2 v 



crit ) 



r b 








1-1/2 





1 



dE = — (1 + ln(2vrA^)) 
2ir 



Thus, the width of the broadest state at the critical point is small in comparison to the total 
length 2N of the spectrum. 
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3.3 Numerical illustration 



In the following, we illustrate the behaviour of the decay widths by results of numerical studies 
for different a and for some finite values of N between 50 and 5000. 



Splitting the sum for finite N in Eq. (35) into its real and imaginary part, one gets 



N 



k-e 



k=-N 



(k — £) 2 + (r/2) 2 



(45) 



and 



A? 



r/2 



1 a £ N (k-sy + (r/2y 



(46) 



The first equation describes the trajectories of the eigenvalues A = £ — |T of TL in the complex 
plane while the second one contains their parametrization with a. We calculated the sum in 



Eq. (45) for different N and fixed values of V and traced their solutions as a function of £. 



In Fig. l.a, the results of the calculations for different N (N = 500,1000,5000) and fixed 
r = 1 are shown as a function of 8. Due to the denominator of the sum, every resonance state 
k with Tfc > 1 has two solutions while there are no solutions when < 1. The number of 
solutions and thus the number of resonance states with > 1 depends on the average slope by 
which the sum approaches the value as a function of £ . 



The result is as follows: For a fixed value of T = 1 there are the more solutions of Eq. (45), 
the larger N. In the limit N — * oo, Eq. ( p5|) can be fulfilled for all resonances, i.e. all resonance 
states have widths larger than an arbitrarily chosen finite value. This confirms the analytical 
result for an infinite number of states, where we have shown that the widths of all states diverge 
at a = a CT it- 



In Fig. l.b, the results of calculations are shown with a fixed number N = 1000 and different 
values of r (r = 0.5, 0.75, 1.0). In this case, the average slope of the different curves is the same 
but the amplitude of the oscillations varies. The larger T the smaller the amplitude is. That 
means, in the case of finite N, the number of resonance states having > T is the smaller the 
larger T. 

As a result, we state the following: When N is a finite number, only a limited number of 
resonance states has widths > T where T is an arbitrarily chosen finite value. Further, the 
larger N, the larger the number of resonance states with > T. On the other hand, the 
larger T, the smaller the number of resonance states with > T. Thus we have two processes 
compensating each other which ensures, that ( |33[ ) holds also in the limit a — > ^ and N — > oo. 

In Fig. 2. a, we illustrate the motion of the eigenvalues in the complex plane for 0.01 < a < 2 
in steps of 0.01 for positive energy £ (the part for negative energy is symmetric to that for posi- 
tive energy). For each resonance state its eigenvalue follows a certain trajectory with increasing 
a. For the lowest values of a, all eigenvalues are near to = k and r^/2 = a. The full line 
gives the estimation for the upper limit of r/2 at a = a cr it according to Eq. (f42|). The estima- 
tion is good in the centre of the spectrum. The deviations at large energies are pure boundary 
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effects. The differences between the different eigenvalue trajectories arise from the finite value 
of N. In the limit of N — > oo all eigenvalues acquire the same behaviour because of the discrete 
translational symmetry on the real energy axis. 

We show in Fig. 2.b the behaviour of all as a function of a for N = 50. At the critical 
value a cr it = 1/vr (indicated by a vertical solid line) the width of the collective resonance state 
k = separates from the widths of the other ones and increases linearly. The slope of To(a)/2 
is equal to 2N + 1 over almost the whole range of a > a cr \t according to Eq. (|32l) . The larger 
slope of To (a) close to a cr i t is a boundary effect and disappears in the limit of N — > oo. 

Fig. 2.c shows Nk° as a function of a for two different values of N. The curves show a sudden 
rise at a cr it and saturate rapidly to 1. The inlet gives a magnification of the curve around the 
critical point which is marked by a vertical solid line. The larger N, the sharper are the changes 
in the slope of Nq(cx). This is a clear numerical indication of the cooperative effect acting over 
the total length of the spectrum. In spite of the fact, that the width Tq of the fast decaying 
state at a ~ a cr it is of the order ln(M), eq. (f44|), its wavefunction carries contributions from 
basis states which are lying far away from the centre of the spectrum and are not overlapped 
by it. These contributions over large energy scales are achieved via the "chain" of overlapping 
neighbouring resonances. As a result, we observe a "macroscopic" order over the whole energy 
scale of the system being much larger than Tq. 

The curve of Nq does, of course, not jump immediately to 100% at the critical point. The 
main reason is the finite number of states taken into account in the calculations. The widths of 
the resonances at the border of the spectrum are in general smaller than those of the resonances 
inside the spectrum. Therefore the chain of neighbouring overlapping states is interrupted at 
energies close to the edges. 

For all the trapped states the corresponding quantity Ni always remains in the order of 
1/M. The wavefunctions of these states are mixed only with those of their next neighbours. 

All the numerical results show that although the phase transition appears mathematically 
for N — ► oo only, the characteristic features of it can be already seen at comparably small values 
of M. 



3.4 Picket-fence level distribution with disturbed translation invariance 

Let us break the translation invariance of the picket-fence model by giving another coupling 
strength to the state in the centre of the system. Suppose: E^ = k und \vy\ = 1 \/k — {0} and 
vq = 1 + D. In this case, the characteristic polynomial (§) reads 

N N N 

Pn(X)= II (fc-A)-ia Yl l[(j-X)-iaDl[(k-\)(k + \) . (47) 

k=-N k=-Nj^k k=l 

Dividing by a convergence ensuring factor and identifying the resulting terms with the product 
representation of sin(a?) and cos(a;), respectively, we get 

sin(7rA) , , . sin(7rA) . . 

P(A) = — - — - + ia cos(ttA) + iaD — = (48) 

7T 7TA 
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in the limiting case N — > oo. In order to study the behaviour of the state in the centre (£ = 0) 
we write A = — ifi and get 

1 



1 

a = — 

TT 



D_ 

7T[l 



+ COth("7T//) 



(49) 



According to this equation, a — > ^ for /x — > oo. The redistribution of the system takes place at 
the same finite value of a cr it = l/vr as in the case of constant coupling. 

We investigate now the behaviour of the system at the critical point, i.e. the type of the 
singularity. Suppose e = 1 — ira and 7T/u = ce~ s with s G R. Using the relation coth(7Tju) — > 1 
for large /j,, we get from ( |49"| ) 

1 = (1 - e)(l + -e +s ) . (50) 
c 

In leading order of the singular part of irn(e), this equation is solved by s = 1 and c = D(l — e) ~ 
D. That means, we have an algebraic singularity: (irfi = ~). The system approaches the singu- 
larity quicker than in the case of a picket fence with translation invariance. 

In the following we like to discuss formula (^) in detail with the help of a numerical il- 
lustration (Figs. 3 and 4). It turned out, that the case with a discriminated state at £ = is 
more complicated. Therefore we have drawn the movement of the complex eigenvalues A in the 
centre of the spectrum for D = —0.5 (Fig. 3. a). We use M = 2iV-|-l = 101 and 301, respectively. 

At a = a cr ;t = l/vr two broad modes arise at the flanks of the spectrum at \£\ ~ 4.5 and 
7.5 for N = 50 and 150, respectively. The larger the number of resonance states is, the larger 
is their distance from the centre. With increasing a the poles of the two resonance states at 
positive and negative energy approach each other in their real part and collide at £ = at a 
certain value a = a c \. At a > a c \, the resonance states remain at £ = 0, one with further 
increasing, the other one with decreasing T. The more resonance states the spectrum contains 
the smaller is a c \, as one can see in Fig. 3.b. Here the imaginary parts of the eigenvalues 
Afc shown in Fig. 3. a are drawn as a function of the coupling strength a. The collision point 
shifts to the critical point of the spectrum (a c \ — ► a cr ;t) if the number M of resonance states 
is enlarged, in spite of the fact that their distance from one another at a c \ is larger when M is 
larger. The two values of a c ± corresponding to M = 101 and M = 301, are indicated by vertical 
dashed lines. The values of a cr it and a C 2 are marked by vertical solid lines. In the limit M — > oo 
the two broad poles appear at £ — > oo. In this limit, a c \ — ► a C rit- At this point of a, the poles 
jump to £ = 0. 

Between a = a c \ and the finite value a = a C 2, there exist three resonance states at £ = 0: 
the two broader poles appearing at a c \ and the original one which is discriminated by the ex- 
ternal coupling by D. The collective mode is one of the two resonance states arising from the 
phase transition at a = a c \. Its imaginary part increases with further increasing a. The other 
broad pole decreases in T with increasing a. Its collision with the discriminated resonance state 
at a = a C 2 shifts both states away from £ = and the imaginary part of both eigenvalues 
decreases. Contrary to the value of a c \ which approaches a cr i t with M — > oo the value of a C 2 
remains almost constant as a function of M. For M — * oo, the value of a C 2 remains larger than 
a cr it- As we will see below, it mainly depends on D. The poles of the trapped states approach 
the values n + 1/2 (with n £ Z) if a — > oo. 
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Figs. 4. a. and 4.b. show the graphs of n \i coth(-/r \i) + D for different D (D = —0.5,0,0.5) 
and fj,/a for several values of a as a function of [i. The points of intersection are the solutions 
of equation ( |49| ) derived under the assumption of an infinite number of states. 

In Fig. 4. a. the coupling parameter is set to a = 0.1 < a cr i t and a = 1/tt = a cr it, respec- 
tively. For the value a < a cr it there exists only one point of intersection with each of the curves 
of coth, lying at small values of \i. Also the state at £ = has a comparably small width in 
the undercritical regime of a and the value of \x increases with increasing D. In the parameter 
range a < a cr j t no broad mode is separated from the other ones. 

At the critical point a = a CT n, where the phase transition takes place, the linear curve fj,/a 
is tangential to the irfx coth(-7r//) + D for D = and is parallel to this function for D = —0.5, 0.5 
(lower and upper thick full line, respectively). So each of these curves has a point of intersection 
with ji/a at = oo. For the case D = —0.5, the intersection at = oo contains two solutions 
for the two broad modes, arising at the borders (at £ = dboo) of the spectrum and colliding at 
£ = 0. Additionally, there is another intersection with \xja at a small value of \x which arises 
from the discriminated state at Eq = 0. 

In Fig. 4.b, we see the same curves 7r//coth(7r//) + D as in Fig. 4. a together with fi/a for 
different values of a > a C rit: « < & C 2, a « a C 2 and a > a C 2- In all cases, the curves have 
intersections at \x = oo. This means, for all values of D a mode exists at £ = with an infinitely 
large width if a > a CT n. The curve 7r/i coth(-7r/i) — 0.5 shows three intersections with \xja in 
the range a cr it < a < a C 2- The two intersections at smaller fi approach each other if a — ► a C 2- 
The value of a C 2 depends obviously on the value of D. In dependence on the negative shift of 
coth(7r//) the slope of the tangent changes and therefore the value of a C 2- At larger values of the 
coupling (a > a C 2) only one solution remains at £ = with fi = oo. 

Finally we have investigated the behaviour of the number Nq of principal components of the 
broad resonance as a function of a for D = ±0.5 and M = 101 states. The results are drawn in 
Fig. 5 together with the former one for D = 0. For all cases the sudden rise of Nq at a ~ a cr it 
can be seen. At the critical point the collective state is created, in all cases, by almost all basis 
states distributed over the whole spectrum. This collective behaviour of a phase transition is 
well pronounced also for D / 0. The curve for D = 0.5 is, however, much smoother than the 
other ones. It rises up not as quickly as the other ones and does not approach the maximum 
value 100% in the range of a shown in the figure. Nevertheless, the characteristics of a phase 
transition can be seen also in this numerical study: Comparing the results for M = 301 states in 
the case of D = 0.5 with those for M = 101 states we see the following tendency. By increasing 
the number of states, Nq for D = 0.5 comes closer to the curves Nq for D = 1 and —0.5. This 
behaviour is a hint to a phase transition (in the next section, Fig. 7. a, we show an example in 
which a phase transition does not occur and the results as a function of an increasing number 
M of states do not show such a tendency). 

3.5 Phase transition 

Summarizing the results of our study with the ideal picket fence, we state the following. 

- For an infinite number of states the imaginary parts of the complex eigenvalues of the 
effective Hamiltonian (||) show a singularity at the finite value a = a cr it = 1/vr. A 
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bifurcation in the widths appears: the width Tj = o of the state in the centre of the spectrum 
increases with further increasing a while the widths of all the other states start to decrease. 
All states but the state at E = are shifted by 1/2 at a = a cr it- For finite M the width 
Tj = o increases linearly with a > a cr it with a slope given by the number M of states included 
in the spectrum. 

- The state in the centre of the spectrum is a collective one in a global sense. The numerical 
studies for finite systems show that it contains components of almost all basic states of 
the system, also of those which are not overlapped by it. This is expressed by the number 
Nf =Q of principal components of the state in the centre of the spectrum which grows, at 
a = Q cr it, suddenly from its minimum value to 100% (maximum mixing). 

- At a = a cr it, the width Tj = o of the state in the centre of the spectrum is much smaller than 
the extension of the spectrum. Therefore, the system does not create locally a collective 
state which traps, with increasing a, further resonance states overlapped by it. 

- At a = a cr ;t, the system suffers a change in its structure: one of the resonance states 
aligns with the decay channel. Its wavefunction collects all the corresponding components 
from the wavefunctions of all the other states (which appear with the same weight in all 
basic wavefunctions because of the symmetry of the problem). Therefore, its width Tj = o 
increases with further increasing a while the widths of all the other states decrease. 

Generally, the behaviour of the order parameter of a system as a function of a control pa- 
rameter characterizes the type of the phase transition. In case of a first order phase transition 
the order parameter shows a jump at a certain finite value of the control parameter a = a cr \t- If 
a higher order phase transition is present, the corresponding derivative of the order parameter 
jumps at a cr it- More precicely, its (n — l)-th derivative jumps in case of an n-th order phase 
transition. 

In our system, the value Tq/M can be understood as the order parameter. It is the width of 
the collective state normalized according to the number of resonance states contributing to its 
formation. This value increases linearly as a function of a with the slope 1/M for a < a CI a and 
with the slope one if a > a cr a. So, the first derivative of Tq/M jumps at a = a CT \t- Therefore we 
conclude, that the formation of a globally collective resonance is a second order phase transition. 
We will prove in the next sections, that this behaviour is universal for all systems showing a 
collective reorganisation. 

In sections 4.5 and 5, we will see that the phase transition is accompanied by an essential 
deviation of the value B (Eq. (|23|)) from 1 in the neighbourhood of a cr it- This means the bi- 
orthogonality of the function system plays an important role in the reordering process. 



4 More realistic systems 

In the previous section, we investigated the properties of the translation invariant picket-fence 
model as a function of the coupling parameter a analytically as well as numerically. The system 
suffers a second-order phase transition at a = ~ which we studied in detail. We consider now the 
behaviour of some more realistic systems as a function of the parameter a in order to answer the 
question whether the results obtained have a general meaning. In detail, energy dependencies 
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and fluctuations within the spectrum will be considered in the following. 



4.1 System with unequally distributed levels and equal coupling strength 

We investigate now the behaviour of systems when the level density is not constant but changes 
as a function of energy. Suppose: = sign(fc) k 2 und = 1 V/c. This means, the level 
density is assumed to decrease linearly with k and to approach the value zero for k — > oo. 

In this case, the characteristic polynomial reads 

N N N 

f N (A)= J] (sign(fc)fc 2 -A)- M ^ [] (sign(j) j 2 - X) (51) 

k=-N k=-N j=-N;j^k 

which can be rewritten as 

Pn(X) =Q N (X)+ia-^Q N (X) (52) 

where Qn{X) = IlfcL-Ar (sign(fc) k 2 — Xj . Dividing Qn{X) by the convergence ensuring factor 
F = -n£Li-(*) 4 weget 



k=l 

In the limit N — > oo, this expression is 



— J_ sin(7rv / A) sinh(7rv / A) = — j sin(-7rv / A) sin(i7rv / A) 



7T^ 7T 

— i 



2tt 2 



cos(7ta/A(1 - i)) - cos(vr\/A(l + i)) (54) 



according to the Weierstrass product representation. In order to find the solution at £ = 0, we 
write A = ifi and get 

} = — cos( 7r \f2jj) — cosh(7rv / 27I) (55) 

F 7T 2 L J 

and finally, with ^Qjv = 

a/S/I cosh(7r^/27I) - cos(7r^/2/7) 

a = — — 7==r ; — 7 i==r- . (ODJ 

7r smh(7r v / 2^J + sin(7r v / 2//J 

The right-hand side of this equation is a monotonically increasing function: a — > oo with — > oo. 
The dilution of the spectrum at large |£| prevents therefore a phase transition. 

Fig. 6. a shows 1^/2 as a function of a for iV = 50 states. A short-lived state is formed, but 
in contrast to the ideal picket fence (Fig. 2.b), no critical value of a can be defined. The width of 
the state in the centre of the spectrum (at £ = 0) increases smoothly as a function of a trapping 
step by step its neighbours. The formation of the short-lived state does not occur by a collective 
interaction of all basis states but by individual trapping of neighboured levels. In other words, 
the short-lived state is not formed by a cooperative effect acting over the whole energy scale of 
the spectrum but is restricted to the energy range overlapped by it. There is no phase transition. 
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The number Nq of principal components in the wavefunction of the state k = is shown 
as a function of a in Fig. 7. a for N = 50 and 150. It supports the conclusion drawn. Also 
this value is a smooth function. There is no hint to a phase transition. The broad state carries 
only components of those basis states which it overlaps. In contrast to the cases where a phase 
transition occurs, the slope of ./Vq decreases with increasing N. The curve ./Vq in Fig. 7. a remains 
smooth unlike the curves for D = 0.5 and M = 101, 301 states in Fig. 5. The collective state is 
created by the basis states of a local energy region overlapped by it. 



4.2 System with unequally distributed levels and unequal coupling strength 

We investigate now the question whether a system with diluted level density at large \£\ shows 
a phase transition if the states at the border are coupled stronger to the decay channel than 
those in the centre of the spectrum. 

To this purpose we rewrite Eq. (M), 



a 



k=l 



+ 



\v-k\ 2 



(57) 



As in the foregoing examples, the spectrum is supposed to be symmetrical (or nearly symmet- 
rical) in relation to E^ = q = 0: — E\. ~ E_^ und Vk ~ V—k- Looking at the solution at E = 0, we 
write A = — Then we get 



from Eq. (57). Since we are interested in the question whether the system shows a phase 
transition at a finite value of a, we have to consider ( |58| ) in the limiting case fj, — > oo and 
N -> oo: 



Ocrit 



lim 

fl— >oo 



N 



\ v k\ 



lim 



(59) 



We approximate Eq. (q9() by an integral, replace the discrete index k by the continuous variable 

x r . Then 



x and assume E\ 



x l und v\ 



= lim 2fx 

"crit ^0 

The integral converges when t 



-Ax 



x l + n 2 
r + l + e \fe > 0. 



lim 2u 2 



--1 



S f + 1 



-ds 



(60) 



Let us consider the following cases. 

1. 2(r + 1) > t. In this case, a cr i t — » with /x — > oo. 

The system is in an overcritical situation for all a > 0. A phase transition does therefore 
not take place. 

2. 2(r + 1) = t. In this case, a = a cr it > remains finite in the limiting case /i — ► oo. 

For a < a cr it, the widths of all states increase. At a = a CT n, a phase transition takes place: 
the short-lived state appears suddenly and a clear separation of time scales, with respect 
to the lifetimes of the resonance states arises, even if the spectrum is infinitely extended. 
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3. 2(r + 1) < t. In this case, a cr it — ► oo with \i — > oo. 

For all finite values of a, there exist states which are not overlapped by the collective 
resonance, whose widths increase with increasing a. Therefore, the state at £ = traps 
new states endlessly. As a consequence, the formation of the short-lived state at £ = 
takes place smoothly. A phase transition does not take place. 

This analytical study shows the following result. To fulfill the conditions for a phase tran- 
sition, the energy dependence of the unperturbed spectra of TL° must be compensated by an 
energy dependent coupling of the individual states to the decay channel. A phase transition 
exists in the cases considered, if the energy dependence of the distribution of the levels Ek is 
opposite to that of the coupling matrix elements v^. According to Eq. (|60|), the critical value of 
a is 

"crit = = 7T~ • (61) 

Otherwise, the system is either in an overcritical regime (corresponding to a cr it — > 0) or in an 
undercritical one (corresponding to a cr i t — > oo). 

Further the width of the broad pole at a ~ a cr it in the compensated case (r ^ 0) can be 
estimated with respect to the length of the spectrum. The characteristic polynomial (Eq. t^j) 
at the energy of the collective state (£ = 0) reads: 

« k h N E k + {T k t N +I*/4 1 ) 

The principal value of that sum gives zero. Therefore we can write: 

i _ r A |Af 

9 2^ 



a 



2 k ^ N (*^) 2 + r=74 

r r N \x\ r , r r Nr+1 <i« 

-ax 



2 J-N (x r+1 ) 2 + T 2 /4 r + 1 Jo s 2 + r 2 /4 



d(4?) 2 /2iV r 

arctan — — — (63) 



r + 1 Jo (^) 2 + l r + 1 V r 

Using the expression (61) for the critical point of the infinite system one gets: 

7T /2A^ +1 \ 

— arctan — — — (64) 



2 v r 

which holds only if T <C 2 • A^ r+1 . In other words, also in the compensated case, the width of 
the fast decaying collective resonance state is much smaller than the extension of the spectrum 
at a coupling strength close to the critical point. 

The compensated case is illustrated in Figs. 6.b and 7.b, the overcompensated one in 6.c and 
7.c. In the compensated case (with r = 1), both the distance between neighbouring levels and 
the coupling strength \vk\ 2 increase linearly with \E\, (Ek = sign(/c) k 2 and \vk\ 2 = \k\ + 1 Vfe). 
In the overcompensated case, the energy dependence of the coupling is chosen stronger than the 
dilution of the spectrum (E^ = sign(/c) k 2 and = \k\ + 1 V7c). Here the coupling increases 
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quadratically whereas the level density decreases linear with E. 



The compensation of the energy dependence of the level density by a corresponding one in 
the coupling strength restores the phase transition (Figs. 6.b and 7.b as compared with 6. a and 
7.a). A collective mode is created by participation of (almost) all basis states. It occurs suddenly 
at a critical value of a. 

In Fig. 7.b Aq for M = 101,301 and 1001 states is drawn. With increasing number M of 
states the curve rises more and more sharply. The critical value is a cr i t = 2/n (indicated by a 
vertical solid line) in accordance with Eq. (|ol"|). 

The maximum value ./Vq = 1 is not reached in Fig. 7.b. Aq is approximately 0.9 and even 
decreases with further increasing a. Drawing the contributions 6oi in the wavefunction <I>o, 



Eq. (21), one sees the following feature. The contributions of the states j coupled more weakly 
to the channel decrease for a > a CT [ t in contrast to those of the states i coupled more strongly: 
Vj < Vi — > \boj\ < \boi\. Since the differences between the Vi are quite large in the cases 
considered, Aq < 1. This holds for both cases, decreasing and increasing energy dependence of 
the level density. In both cases, the coefficients \boj\ are spread at large a, e.g. at a « 4a cr it- 
This is in contrast to the case of the ideal picket fence with equal coupling strengths, in which 
all coefficients approach the value 1/M for a > 2a cr it. 

In the case of an overcompensation (Fig. 6.c and 7.c) the critical point is shifted to very 
small values in accordance with a cr it — > 0, Eq. (|60[). The number of principal components of the 
broad mode jumps up to 75%. Then it decreases and saturates at around 57%. 

In an additional calculation, we bounded the spectrum from below: we investigated the case 
with Ek = k 2 , \vk\ 2 = k + lV/c, k > 0. Also in this system, a phase transition takes place at 
a cr it = 2 /it as in the case shown in Figs. 6.b and 7.b. In all cases, the broad mode appears in 
the energetical centre of the spectrum. 



4.3 System with unfolded Gaussian distributed levels 

It is interesting to learn whether the conditions for a phase transition must be fulfilled strictly 
or only on the average. In order to answer this question in the affirmative, we perform the 
following numerical analysis. We choose an unfolded Gaussian-Orthogonal Ensemble (GOE) for 
the distribution of the eigenvalues of TL° and a Gaussian distributed coupling vector V with 
mean value (v) = 1 and variance Av = 0.01. 

The decay widths as a function of the coupling parameter a are drawn in Fig. 6.d. for 
A = 50. The broad mode separates from the other ones at approximately a CT a = 1/ir with a 
slope of 2 A + 1. The features of the phase transition are not as clearly pronounced in this figure 
as in the case of the ideal picket fence. Nevertheless, the differences to Fig. 6. a., where no phase 
transition occurs, are obvious. Even for the comparably small number of states (M = 101), the 
fluctuations in the distribution of the levels and the coupling vector do not destroy the nature 
of the reorganisation process. 

As shown in the foregoing sections, the number of principal components A^ is a sensitive 
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quantity to measure the global collectivity of the separation process. In Fig. 7.d. Nq of the 
collective mode at E = is drawn as a function of a for N = 50, 150, 250, 500. For increasing 
M = 2N + 1 the curves rise up more suddenly and the slope near a = 1/n gets steeper. All the 
curves approach the maximum value of ./Vq very fast for values a > a cr \f 

The features of the second-order phase transition are better expressed if more resonance 
states are considered. For large N the irregularities in the distribution of the Ej- and Vk are 
almost unimportant. This proves, that the conditions derived in the former sections have to 
be fullfilled only on the average. Also in the ergodic case of a GOE distributed spectrum, the 
reorganisation of the spectrum can be understood as a second-order phase transition. The fluc- 
tuations within the spectrum will be washed out if the conditions for a phase transition derived 
from Eq. ([k]) are fulfilled on the average. 



4.4 System with complex coupling parameter a 

Up to now, we considered the system to be described by the Hamiltonian TL, Eq. (|5|), where 
the coupling between system and continuum is supposed to be real and TL is non-hermitian. 
There may be an additional part [5 VV + in the coupling term by which a collective state of 
another (internal) type is created. This collective state is shifted by an energy AS from the 
group of the remaining A — 1 states @. The structure of both parts VV + and VV + is the 
same. The difference is the non-hermiticity of the external coupling term in the first case and 
the hermiticity in the second case. 

We are interested in the question whether the additional term has an influence on the phase 
transition. Investigating this question, we restrict ourselves to the case VV + = VV + , i.e. an 
angle zero between the vectors V and V. Further, the characteristic polynomial (JgJ) does not 
contain, in the one-channel case, the phases of the coupling matrix elements, but only \vk\ 2 - It 
is justified, therefore, to replace a by a + if3 in the equations considered in the previous sections 
in order to get conclusions on the influence of the term with (3 7^ on the phase transition. 

Considering the picket-fence model with equal coupling strength (which we studied in section 
3 for P = 0) Eq. fl26|) must be replaced by 




P(A) = sin(vrA) + iir(a + iff) cos(vrA) (65) 
Using the representation A = £ — iV, one gets 

-an „7rr\ \ /„„„/„z?\\ 

(66) 

This equation has a solution, when the determinant of the matrix vanishes. This condition gives 

1 ( ^g + lf + ^p)^ 
T = ^ ln {(na-iy + (npr) " (6?) 

Eq. ( |67| ) has no singularity when P ^ 0. This means, the singularity in the widths of the reso- 



nances, obtained for (3 = in Eqs. (27) vanishes when P ^ 0. 
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This result can be understood as follows. In [^] the distribution of exceptional points is 
investigated for systems described by a Hamiltonian of the type H = H° + 5cH . The excep- 
tional points of such a system are those points in the parameter space of a at which two (or 
more) eigenvalues coincide (for a more detailed discussion see for example Jl6|]). The coupling 
constant a can be a real, imaginary or, more generally, a complex number. The distribution 
of the exceptional points is determined by the matrices H° and H 1 . It is independent of the 
value of a which determines, for its part, the positions of the (in general complex) eigenvalues 
of H. In the one-channel case (in which H 1 has rank one), there exist M — 1 exceptional points 
corresponding to the crossing of the collective state with each of the other M — 1 states. 

For systems which show a phase transition, allM — 1 exceptional points converge to the finite 
purely real value of a = a cr \t in the limit M — > oo. For finite systems, almost all exceptional 
points are near to this accumulation point. This result holds not only for the ideal picket fence, 
but for all systems, which suffer a phase transition [15]. 



For the systems investigated in the present paper it follows. The accumulation point being 
determined by the matrices H° and VV T , is independent of the coupling parameter a. When 
the system with a purely imaginary coupling (i.e. (3 = 0) shows a phase transition, all (M — 1) 
exceptional points are met if a approaches the critical value. The collective mode repels with 
all the other ones simultaniously, i.e. all states run through their exceptional point at a = a cr \t- 
This means all M — 1 exceptional points are accumulated at a cr it- In that case, (<&j|<&j) diverges 
for all i simultaneously in the limes a — > a cr it- In fact, the dimension of the eigenspace collapses 
from N to 1 = \&j) Vi, j) if a hits the accumulation point. Therefore, also B diverges at 

a cr i t . If the coupling parameter is complex (/? / 0), however, the system passes the accumulation 
point in a certain distance in the complex (a, /?)-plane. As a result, the singularity at a cr it will 
be avoided. For |a| ~ |a C rit|) the quantity B does not diverge if M — > oo but reaches a certain 
maximum value. Refering to this result, we claim that, according to a rigorous mathematical 
definition, the phase transition will be destroyed by any given nonvanishing real part (3 in the 
coupling parameter (a detailed discussion of this aspect is given in fll5|| ). 



Let us illustrate this result by means of a numerical study. To that purpose, we replace 
a by a ■ e lip . We choose M = 2N + 1 = 101, as usually, and perform the calculations for 
ip = 1°, 10°, 45°, 80°, 89° by varying a. The eigenvalues of Tt° (for ,9 = 0) and the coupling 
matrix elements are chosen to be Ek = h and Vk = 1, respectively. In Fig. 8. a. we have drawn 
the number of principal components of the collective resonance state as a function of a. 

The numerical results show a clear difference between the cases with small and large angle 
ip. The larger (p, the less is the number of basis states contributing to the collective state at a 
given a > a C rit- Further, the curves rise up more smoothly when ip is larger. For large angles 
the maximum value A^q = 1 is not reached at the maximum value a = 2 shown in the figure. 
Thus, in the case of the finite spectrum studied numerically, the reorganisation process is getting 
smoother the larger (p. In other words, the reorganisation process is washed out, if the system 
cannot hit the accumulation point of the exceptional points, but has to pass it in a certain 
distance in the complex parameter space. 

This conclusion is underlined by the results given in Fig 8.b. Here we have fixed the angle ip 
to 45° and varied the number of states included in the spectrum (M = 101 and 1001 states). As 
one can see, the characteristic features of the curves are not changed by changing the number 
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of basis states. The transition remains smooth also for N = 1001 states over the whole range 
of a. So, these systems can not be characterized by a critical point. We see a critical region of 
a, which will be the larger, the larger the nearest distance between the accumulation point and 
the eigenvalues in the complex plane is. 



5 Discussion of the results 

The analytical and numerical investigations represented in the foregoing sections point to simi- 
larities and differences in the behaviour of the different systems under the influence of varying 
strength of the coupling to the continuum (decay channel). In any case, a restructuring in the 
system takes place (or starts to take place) when the coupling parameter a is large enough. A 
collective state which is aligned with the decay channel is formed in the centre of the spectrum. 
Its wavefunction is coherently mixed in the set of basis wavefunctions of the corresponding closed 
system. The trapped states have incoherently mixed wavefunctions. Beyond a certain value of 
the control parameter a, two different time scales exist (bifurcation of the widths). 

In some cases, the restructuring in the system can be identified as a second-order phase 
transition. The separation of different time scales occurs suddenly at a critical value a cr i t and 
is a collective effect of the whole spectrum. The order parameter Tq/M increases linearily as a 
function of the control parameter a with a universal slope one, as soon as the control parameter 
is larger than its critical value. In other cases, the separation of time scales occurs successively by 
individual trapping of neighbouring resonance states. In that case, the collectivity is restricted 
by the extension of the energy region, overlapped by the fast decaying resonance state. The 
process of reorganisation continues up to a — > oo. 

The differences in the behaviour of systems which show a phase transition to those which 
do not, can be seen nicely in the example of the ideal picket fence. In the case with equally dis- 
tributed levels coupled with the same strength to one common channel, all states are equivalent. 
Consequently, the direction of the energy shift accompanying the local resonance trapping is 
undefined and the local resonance trapping is hindered. The redistribution of the system under 
the influence of the decay channel can take place only collectively. The quantity Tq/M rises 
linearly in a > a cr j t with slope 1. 

More realistic systems show a phase transition when the energy dependence of the level den- 
sity is compensated by an energy dependence of the coupling strength. For example, a dilution of 
the level density can be compensated by a corresponding enhancement of the coupling strength. 
The critical value a cr i t is well determined. Further when the system is bounded from below, a 
phase transition occurs under the same conditions as for non-bounded systems. It occurs at the 
same critical value (r + I) /it. 

Moreover we could show, that the conditions for a phase transition do not have to be fulfilled 
strictly, but only on the average. Small irregularities in the energy dependence of the levels or in 
the distribution of the coupling vectors are washed out, if the number of states in the spectrum 
is sufficiently high. Also the generic case of GOE-distributed states shows the features of a phase 
transition even for a comparably small number of states. The fast decaying state is created by 
all states of the spectrum, independently of whether they are overlapped by it or not. 
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Another characteristic feature of the phase transition is the mixing of the wavefunctions. In 
the case of the ideal picket fence, it changes suddenly at a cr it from its minimum value Nq = 1/M 
to the maximum value A^q = 1 for the state i = 0. The width of the collective resonance state 
at a = Q cr it is of the order of ln(M) whereas the extension of the spectrum is equal to M. Also 
in the more general case, where the energy dependence of the level density is compensated by 
the coupling strength of the resonance states, we could prove, that the width of the collective 
state is much smaller than the extension of the spectrum for couplings close to the critical point. 
Nevertheless, the collective state carries contributions of (almost) all basis states even if they 
are, in the case of a finite spectrum, close to the borders. When a ^ a C rit> the value of Nq is 
independent of a. This holds also if the system is disturbed by random pertubations, where the 
compensation conditions are fulfiled only on the average. 

The situation is different if the energy dependence of the level density is not compensated 
by the energy dependence of the coupling vector. In such a case, the local resonance trapping 
between neighbouring states is not hindered but occurs successively starting in the region with 
the largest level density or coupling strength. If the level density is larger in the centre of the 
spectrum than at other energies and the coupling strength to the channel is the same for all 
levels, the collective mode is created only locally for any finite value of a. This local collective 
state traps successively more and more resonance states in direction to the border of the spec- 
trum. This process of local resonance trapping with increasing width of one state in the centre 
continues endlessly up to a — > oo (in the limit M — >• oo). Although the short-lived state has 
collective properties it is not created by all basis states of the spectrum but only by those which 
are overlapped by it. There is no phase transition at all. 

The resonance structure of the system is, in the cases considered, symmetrical in relation 
to the critical value a cr it of the control parameter although the number of long-lived resonance 
states for a < a cr j t and for a > a cr it differs by 1. As an example, the picket-fence distribution 
with level distance 1 and equal coupling strength to the continuum remains a picket-fence dis- 
tribution also at a > a cr it but is shifted in energy by 1/2 of the level distance. 

For finite M, collective states may be caused also by an additional real part to the Hamilto- 
nian, e.g. 7i' = TC° + (3VV + , leaving Ji! hermitian. It is called internal collectivity in contrast 
to the external collectivity discussed above. In such a case, the process of reorganization in the 
system occurs smoothly. A phase transition does not take place. The eigenfunctions of Ji! are 
orthogonal in the usual manner: (<&£|<I>^) = 1 for all i and a. 

Also the non-hermiticity of the Hamiltonian H is, however, not sufficient for the appearance 
of a phase transition, as the results presented in the foregoing sections show. In any case, the 
value B characterizing the bi-orthogonality of the set of eigenfunctions of TL plays a decisive role. 
Only when it becomes essentially, i.e. when B 3> 1 (see Eq. (p3|)) at a certain well-defined value 
of a, a phase transition takes place. When, however, the reordering of the system takes place suc- 
cessively in a limited region of the spectrum with n <C M states, then B^ = ^ J27=ii^i\^i) > 1 
but B is close to 1. In this case, the reorganisation in the system does not occur collectively but 
smoothly as a function of a. 

As a result, in the case of a phase transition the bi-orthogonality of all the eigenfunctions 
of TC is maximal at (almost) the same value of a and, according to Eq. (p3[), B 3> 1 at a C rit- 
For illustration we show in Fig. 9 the value of B as a function of a for four different cases. 
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The theoretical value of a cr it is marked by a vertical solid line. Only in the case without phase 
transition, this sum is always very close to 1 while it has a clearly expressed maximum at the 
critical point a CT it whenever a phase transition occurs. Further, the eigenfunctions of 7i are 
orthogonal in the usual manner for a <C a C rit as well as for a 3> a C rit- Here, B ~ 1. 

This result can nicely be illustrated by means of the exceptional points defined as the cross- 
ing points of resonance states. They are determined by the structure of the different parts of 
the effective Hamiltonian (see sect. 4.4) but are independent of the coupling parameter. In the 
one-channel case considered by us, they accumulate at one point in the limit M — > oo. In the 
case the system shows a phase transition, the eigenvalues meet this accumulation point when 
considered as a function of the coupling parameter and B — > oo in its neighbourhood. 

Therefore, we may differentiate between four situations: (i) the exceptional points accu- 
mulate at a finite real value in the complex (a, /3)-plane and the system goes through the ac- 
cumulation point, (ii) the exceptional points accumulate at a finite real value in the complex 
(a, /3)-plane but the phase (p of the coupling hinders the system to hit the accumulation point, 
(iii) the exceptional points accumulate at a, (3 = 0, (iv) the exceptional points do not accumulate 
at all but they are spread over the whole complex (a, /3)-plane with diverging absolute value of 
the coupling parameter. Examples are (i) the compensated case, (ii) the system with complex 



coupling, (iii) the overcompensated case, (iv) the undercompensated case (for details see [15]) 



The stochastic processes described by the (complex) partial widths are much larger for 
a a cr j t than at other values of a. This is expressed by the relation Tj = ji c / where ji c 
is the partial width of the state i in relation to the (only) decay channel c Since ^> 1 

near a cr it, it follows r, <C ji c for a a C rit- The structure observed in the cross section is deter- 
mined by the S-matrix, Eq. (Q). It depends essentially, according to Eq. (Q), on the length of 
the spectrum (i.e on the values 7j c of all the states), but not on the width Tj = o of the collective 
state. For illustration, |1 — Su\ 2 is shown in Fig. 10 for three different values of a > a C rit- The 
width of the collective state is To/2 = 0.84, 1.09 and 19.9, respectively, for the three values of 
a considered. The width Tq has almost nothing in common with the structure observed in the 
cross section as one easily sees from the figure. 

The same result follows also from our analytical considerations. At the critical point, the 
sum of the widths of all states is smaller than the total length M = 2N + 1 of the spectrum by 
a factor tt according to Eq. (|33|). Furthermore, the width of the broadest state (in the centre of 
the spectrum) is in the order of In N in the ideal picket-fence model. In the limit of large N it 
is in general tiny compared to the length of the spectrum according to Eq. (|64|). This shows, 
that the transition is caused by the cooperative behaviour of all states. It is not caused by the 
overlap of the complete spectrum by one of the states. In the cross section, we see a structure 
of the extension of the length of the spectrum since all states are coupled to the decay channel. 
The width of the broadest state is much smaller than this structure. 



The numerical results show further that the number M of states need not necessarily 
be infinite in deciding the question whether the transition is of second order or not. The 
second-order phase transition is well expressed already for a relatively small number of states 
(M = 2N + 1 = 101 up to 1001 in our calculations) in all cases in which the analytical study 
shows a phase transition in the limit M — > oo. 
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In our analytical and numerical studies the limiting case M — > oo is achieved by an extension 
of the length of the spectrum. It is worthwhile to note that the results are the same if, instead, 
the length of the spectrum is kept fixed at some finite value and the level density approaches oo 
with M — > oo. 

In many-particle systems, the level density depends on energy. In nuclei, it increases expo- 
nentially with energy. The coupling strength of the states to the continuum decreases, however, 
with energy due to the increasing contribution of many-particle many-hole configurations to the 
wavefunctions of the states ("compound nucleus states"). It is an interesting question whether 
in such a system the increasing level density is " compensated" in a certain energy range of the 
spectrum by the decreasing mean coupling strength so that the condition for a second-order 
phase transition is fulfilled. 
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Figure captions 

Fig. 1. 

Numerical illustration of Eq. (^) for different values of A and T in the positive energy range, (a) 
r = 1,N = 500 (full line), 1000 (dashed line), 5000 (dash-dotted line), (b) A = 1000, T = 0.5 
(full line), 0.75 (dashed line), 1.0 (dash-dotted line). For details see text. 

Fig. 2. 

(a) The motion of the eigenvalues Xk in the complex plane with increasing a for the ideal picket 
fence (N = 50). The ordinate represents the values of r^/2 while the abcisse those of E^. The 
full line shows the behaviour of Eq. (f4l|). Only a part of the spectrum is shown, (b) Tfc/2 as a 
function of a for A = 50. (c) Nq as a function of a for N = 50 (full line), 150 (dotted line). 

Fig. 3 

(a) Eigenvalues = £/. — i/2 in the complex plane for a small energy range around the 
centre. Ek = kMk and Vk = 1 Vfc but vq = 0.5. N = 50 (rhombs), 150 (crosses), (b) r& as a 
function of (a). = kMk and = 1 \/k but vq = 0.5. N = 50, 150. 

Fig. 4 

Numerical illustration of Eq. (f49|). 7r/i coth(7r/i) + D (left ordinate scale) for D = —0.5,0,0.5 (full 
lines) and [ija (right ordinate scale) for different a. (a) for a = 0.1 (dashed line), 1/tt (solid 
thin line), (b) for a cr i t < a < a C 2 (solid thin line), a ~ a C 2 (dashed line), a > a C 2 (dash-dotted 
line). For details see text. 

Fig. 5 

Number Aq of principal components as a function of a for the three different picket-fence dis- 
tributions with vq = 0.5, A = 50 (dashed line), vq = 1, N = 50 (full line), and vq = 2, A = 50 
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(dash-dotted line), 150 (double-dotted-dashed line). 
Fig. 6 

as a function of a for N = 50. (a) Ek = sign(fc) k 2 , v\~ = 1 Vfe, (b) = sign(fc) A: 2 , v 2 = 
|fc| + 1, (c) = sign(fe) A; 2 , Vk = \k\ + 1, (d) E^ distributed according to an unfolded GOE, Vk 
are Gaussian distributed with mean value 1 and variance 0.1. 

Fig. 7 

Nq as a function of a for (a) Ek = sign(A;) k 2 , = 1 Vfc, iV = 50 (full line), 150 (dashed line); 
(b) £ fc = sign(A;) /c 2 , v\ = \k\ + 1, N = 50 (full line), 150 (dotted line), 500 (dashed line); (c) 
Ek = sign(/c) k 2 , = \k\ + 1, N = 50 (full line); (d) E^ distributed according to an unfolded 
GOE, Vk are Gaussian distributed with mean value 1 and variance 0.1. N = 50 (full line), 150 
(dotted line), 250 (dashed line), 500 (dash-dotted line). 

Fig. 8 

A^q as a function of a for a system with complex coupling. Simultaneously, (3 changes accord- 
ing to (3 = atan(tp). (a) The values ip = 1°, 10°, 45°, 80°, 89° are shown as full, long-dashed, 
short-dashed, dotted, dash-dotted line, respectively. N = 50. (b) ip = 45°, iV = 50 (full line), 
N = 150 (dotted line). 

Fig. 9 

The value of B, Eq. (|23|), as a function of a. The thick full line shows the case Ef. = 
sign(fc) k 2 , Vk = 1 Vk,N = 50, thin solid line the ideal picket fence for N = 50, dashed line 
the ideal picket fence for ./V = 150, dash-dotted line the unfolded GOE as in Fig. 7.d for iV = 150. 

Fig. 10 



|1 - Sii| 2 for a = 1/vr (full line), l/vr + 0.05 (dashed line), 1/vr + 0.1 (dash-dotted line). For 
comparison: r /2 = 0.84 for a = l/vr, r /2 = 1.09 for a = 1/vr + 0.05 and r /2 = 19.9 for 
a = l/vr + 0.1. 
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